Detecting very-high-frequency rehc gravitational waves 
by electromagnetic wave polarizations in a waveguide * 

■ Ming-Lei Tong and Yang Zhang^^ 

I Centre for Astrophysics, University of Science and Technology of China, Hefei 230026 

(N 

> 
O 

o 

en I The polarization vector (PV) of an electromagnetic wave (EW) will experience a 

rotation in a region of spacetime perturbed by gravitational waves (GWs). Based on 
qh' this idea. Cruise's group has built an annular waveguide to detect GWs. We give 

detailed calculations of the rotations of the polarization vector of an EW caused by 
incident GWs from various directions and in various polarization states, and then 
analyze the accumulative effects on the polarization vector when the EW passes n 
^ ■ cycles along the annular waveguide. We reexamine the feasibility and limitation of 

■ this method to detect GWs of high frequency around 100 MHz, in particular, the relic 
Q^ ■ gravitational waves (RGWs). By comparing the spectrum of RGWs in the accelerating 
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universe with the detector sensitivity of the current waveguide, it is found that the 
amplitude of the RGWs is too low to be detected by the waveguide detectors currently 



. running. Possible ways of improvements on detection are discussed also. 



' 1 Introduction 



GW is one of important predictions of general relativity. Although there has been 
an indirect evidence of GW radiation from the binary pulsar B1913+16 [1], so far 
direct detection of GWs have not been accomplished yet. GWs can have different 
frequencies generated by various kinds of sources. Currently, besides the conventional 
method of cryogenic resonant bar [2] , a number of detectors using new techniques have 
been running or under construction aiming at direct signals of GWs. For a frequency 
range 1 ~ 10^ Hz, the method of ground-based laser interferometers applies, such as 
LIGO [3], Virgo [4], and TAMA [5], . For a lower frequency range 10"^ ~ 1 Hz, the 
space-based laser interferometers can be used, such as LISA [6] under planning. For 
much lower frequencies ~ 10"^^ Hz, detections of CMB polarization of "magnetic" 
type would also give direct evidence of GWs [7]. There have also been attempts to 
detect GWs of very high frequencies from MHz to GHz, employing various techniques, 
such as laser beam [8]. One interesting method proposed by Cruise uses linearly 
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polarized EWs [9, 10]. When GWs pass through the region of waveguide, the direction 
of PV of EWs will generally experience a rotation [9] . A prototype gravitational waves 
detector has been built by Cruise's group [10], which consists mainly of one, or several, 
annular waveguide of a shape of torus. As a merit of this method, depending upon 
the size of the waveguide, GWs in a very high frequency range 10^ ~ 10^ Hz can be 
detected, which is not covered by the laser interferometer method. Note that the GWs 
in the frequency range 10^ ~ 10^ Hz are generally not generated by usual astrophysical 
processes, such as binary neutron stars, binary black holes, merging of neutron stars 
or black holes, and collapse of stars [11] [12]. However, the background of RGWs has 
a spectrum stretching over a whole range of 10~^^ ~ 10^^ Hz [13, 14]. Depending 
on the frequency ranges, its different portions can be detected by different method. 
For instance, the very low frequency range 10^^^ ^ 10^^^ Hz can be detected by the 
curl type of polarization in CMB [7], the low frequency range 10~^ ~ 10~^ Hz can be 
detected by LISA, the mediate frequency range 10^ ~ 10^ Hz is covered by LIGO, and 
the very high frequency range 10^ ~ 10^ Hz can be the detection object of Cruise's 
EWs polarization method. Therefore, one of the main object of detection by the 
annular waveguide is the very high frequency RGWs. The detection of high frequency 
RGWs from MHz to GHz is in complimentary to the usual detectors working in the 
range of 10^'^ ~ 10"^ Hz. RGWs is a stochastic background that are generated by 
the inflationary expansion of the early Universe [13, 14, 15, 16, 17], and its spectrum 
depends sensitively on the inflationary and the subsequent reheating stages. Besides, 
the currently accelerating expansion also affects both the shape and the amplitude of 
the RGW spectrum [13, 14, 17]. RGWs carry take a valuable information about the 
Universe, therefore, their detection is much desired and will provide a new window of 
astronomy. 

In this paper we give a comprehensive study of the rotations of PV of EWs in a 
conducting torus caused by incident GWs, and explore the feasibility and limitation 
of Cruise's method of detecting GWs by polarized EWs in the annular waveguide. 
Firstly, we briefly review the RGWs in the currently accelerating universe. Secondly, 
we shall present detailed calculations of rotations of the PV of EWs in the waveguide 
caused by the incoming GWs from from various directions and in various polarization 
states, thereby we analyze the multiple-cycling accumulating effect and the resonance 
when the circling frequency of EWs is nearly equal to that of GWs. Thirdly, we shall 
examine the possible detection of the RGWs by the annular waveguide system around 
100 MHz, comparing the predicted spectrum of RGWs in the accelerating Universe 
with the sensitivity of the detector [10]. Finally, we give the conclusions and possible 
ways of improvements for detection. 

2 Relic gravitational waves 

In an expanding universe RGWs can be regarded as small perturbations to the Robertson- 
Walker metric, 

ds^ = a^{T)[-dT'^ + {dij + hij) dx'dx^], (1) 



where a(r) is the scale factor, r is the conformal time, and hij is transverse-traceless 

di = 0, 8'^ hij = 0, (2) 

representing RGWs. Among six components h^j there are only two independent (two 
polarization states). Generally, \hij\ <C 1. The wave equation for RGWs is 

df,{V^d^hij{^,T)) = o. (3) 

The solution hij of Eq.(3) and the spectrum h{v,TH) defined via 

<h'^hij>= / h\k,TH)-r- (4) 



have been given for an accelerating universe [14, 17, 18]. Fig.l plots r/f), which 
depends on the accelerating parameter 7, the inflation parameter /?, the reheating 
parameter f3s, the tensor /scalar ratio r, and the redshift ze sX the time te of equality 
of dark energy and matter is given by 

l + ZE = ^^^{-^)^. (5) 

Since the annular waveguide is to detect RGWs of frequencies ~ 10* Hz, we quote the 
analytic approximate spectrum in this range [14]: 

where k is the comoving wavenumber related to the physical frequency by = 2TTa(TH) ' 
^0 is a constant determined by the CMB anisotropics [7, 14, 18, 19], e = (1+/3)(1 — 7) /j 
is a small parameter, and kn = 2tt^ and kg — W^^kn- Note that the RGWs hij 
described above exist everywhere and all the time in the Universe. We may simply say 
that the Universe is filled with a stochastic background, consisting of all the modes 
of different wave-vector k'^={k'^, k^, k'^, k^). So the RGWs serve as an object for GWs 
detections. 

In the frequency range ~ 100 MHz for the waveguide detector, RGWs can be 
approximated as plane waves. A beam of monochromatic plane GWs with a wave- 
vector can be generally written as the following form [20] 

hij = Re{Aije'% (7) 

where Aij represents the amplitude and 4> is the phase of GWs, 

^ = kf,x^ = gi^.k^x", (8) 

with x^^ being the point of spacetime that the waves pass. 



3 The Annular Waveguide 

Consider an annular waveguide of shape of a torus, as shown Fig.2. Its radius is R, 
and the cross section is a rectangle with sides a > b, both being much less than R, 



with parameters 
r = 0.22 
nA = 0.75 
Ps = -0.3 
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Figure 1: The spectrum h{v,r]H) of RGW in the accelerating universe. 



say, o, 6 ~ 1 cm, and i? ~ 1 m. (Note that the waveguide actuahy employed by 
Cruise's group [10] actuahy has a shape of rectangle, instead of a torus. For simplicity 
of analysis, here we consider a torus since the working mechanism is the same.) Inside 
the torus, one can input a beam of linearly polarized EW propagating around the 
toroidal loop, and the beam consists of a TEio mode (transverse electric field) with 
the electric held E pointing along the x'^— axis. The EWs are a microwave, say, with 
a wavelength Ag 1 mm and a frequency Ve = cjXe ^ 10^^ Hz. The guided EWs of 
TEio mode in the torus travel around the loop at a group speed [21], 



(9) 



where c is the speed of light, and v is very close to c. For instance, for Ag ~ 1mm 
and a = 1cm, the difference between the two velocities is ~ {\e/2af/2 ~ 10^3. As 
is well-known, for a TEio mode to exist in the waveguide, one has to Ae < 2a. The 
angular velocity of EWs cycling around the loop is then 




(10) 



As will be seen later, when the angular frequency lo of the incident GWs is very close 
to the i-6-, at the resonant condition, the detector responds most sensitively to 
the GWs. Therefore, such a device of given radius R will primarily detect GWs of a 
resonant frequency around 

(11) 



2'kR 



For example, if the radius is i? = 1 m, then the frequency of GWs to be detected is 
t'g ~ 5 X lO'^ Hz, some two orders smaller than fg. As a merit, by adjusting the size R, 
the frequency of GWs to be detected can vary accordingly. For the 4-dim spacetime, 
one can choose a coordinate system {x^} with /u = 0, 1, 2, 3 and = ct, such that the 
waveguide lies on the {x^^x^) plane as shown in Fig. 2. Note that the geometric size 
R of waveguide is negligibly small in comparison with the Hubble's radius ~ c/H, so 
that the effect of cosmic expansion on the torus can be totally neglected. 



A 



The waveguide 




Figure 2: A sketch map of the annular waveguide. The cross section of the waveguide is 
rectangle with a and h being much less than the radius R. EWs travel inside the waveguide 
and GWs propagate along the x^-axis. 

Let us consider a beam of GWs passing through the detector. Assume that the 
wavelength Ag of the GWs is much longer than the wavelength Ag of the EWs in the 
waveguide, i.e., Xg » Ag, so that the geometric optics approximation applies in de- 
scribing the EWs [9, 20]. In fact, this assumption on the incident GWs is automatically 
satisfied if the GWs satisfy the resonant condition. The PV of the linearly polarized 
EWs can be described by a 4-vector, IT^ = (II'^, H^, H^, H^), which is real and normal 
to the wave vector P^^ of the EWs 

n^P" = 0, (12) 

and satisfies the normalized condition [20, 9] 

n^n^ = 1. (13) 

Eq.(12) tells that one can add a multiple of to 11'^ without affecting any physical 
measurements [20], since P^ is a null vector with P^P'^^ = 0. Suppose that the EWs 
are propagating along the x^— axis with the wave vector P^ = (P'', P^, 0, 0), which 
satisfies P^P'^ = 0. Then, by Eq.(12), the PV of EWs can be generally written as 
jj/. ^ (KpO^Kpi^n2,n3), where K is an arbitrary constant. Then, Eq.(13) leads to 

|n2|2 + |n3|2 = 1. (14) 

Since initially the electric field E of the EWs inside the torus is set to be along the 
x^— axis and 11* is, by definition, in the direction of E, so the initial PV is 

n'^ = (0,0,0,1), (15) 

i.e., initially the PV has a vanishing 11^ = component. 

However, the presence of GWs will cause a rotation of the IT^ about the direction 
of propagation, generating a non- vanishing 11^ 7^ i.e., a component E'^ 7^ of the 



electric field E of the EWs in the waveguide. One puts an electric field probe inside the 
waveguide at the intersection of the x^— axis and the torus. The probe is on the line of 

axis, so that it can probe the non-vanishing electric field i?^ due to the rotation of 
E caused by the GWs [10]. The GWs induce an electric voltage on the electric probe 
V = E^al sin (27rfei), where Eq is the TEio mode electric field in the waveguide, / is 
the length of the conducting probe. 

In the geometric optics approximation, the motion of 11'* is described as being 
parallel-propagating along the rays of Ews with the equation 

— ^Ti,U^—^0, (16) 

where s is an affine parameter, which can be chosen to he s = t/To with Tq = l-KjujQ 
being the period of the EWs travelling around the torus. Note that when s goes from 
to 1, the EWs go one cycle around the torus. Since a, 6 <C i?, one can view the EWs 
in the waveguide as travelling along the 1— dimensional loop path 

, - sin 27rs, cos 27rs, j , (17) 

where v is the group speed of the EWs. 

With the initial setup of the polarization of EWs in the torus, we only need to 
consider the component of the polarization in the following. 



4 Change of II^ 

Even the setup of the waveguide detector is fixed in laboratory, GWs propagating 
in space may come in any direction randomly. Therefore, we need to determine the 
rotation of polarization of EWs caused by GWs travelling along the directions x*, 
z = 1,2, 3, respectively. 

A. GWs travelling along the positive x^-axis 

Consider a beam of monochromatic plane GWs travelling along the positive x^- 
axis with a wave vector = (27r/Ag, 0, 0, 27r/Ag). As the GWs just pass the annular 
waveguide whose position is given by Eq.(17), then substituting it into Eq.(8) yields 
the phase of the GWs at the point inside the annular waveguide 

(f) = —2ttsuj/uo. (18) 

Here u = l-KcjXg represents the angular frequency of GWs, and = {^-kITq) is the 
cycling angular frequency of EWs around the torus. As mentioned before, one can 
take the flat spacetime slightly perturbed by GWs to represent the local region of the 
waveguide. In the transverse traceless (TT) gauge, the metric tensor can be written 
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and 



where /i® and hi^ denote the + and x modes of polarization of GWs, respectively. In 
general, these two modes of GWs may be not coherent, i.e. their phases are random 
and independent, a situation similar to the natural light of EWs. If the + and x 
modes have the same phase 4>, which is called as the linearly polarized GWs [20] , then 
by Eq.(7) one has 

/i© = cos 4>, = cos (f), (19) 

where ^0 and ^ig, are real numbers. This case of a linearly polarized GWs will be 
discussed in the following, otherwise we give a clear indication. 

As can be checked, the change in 11"^ due to GWs is of order of so in the 

subsequent calculation, 11^ = 1 is assumed. To calculate the change of 11^ up to the 
linear order of hij, one needs the following Christoffel components 

rii = -ttA^ sin (t)/Xg, 

ri2 = 7rAesin(^/Ag, (20) 

other components are either zero or of order having no contributions. Integrating 

Eq.(16) gives the expression of the change in 11^ around one circle of the torus 

Substituting Eqs.(17) and (20) into the integration, one has 

An^ = J (^A^ sin ^27rs— ^ cos 27rs - A^ sin ^27rs— ^ sin27rs^ ds. (22) 

Carrying out integration yields the result 



Atf = ^ (1 - cos {2^^)) - ^ sin {27rw)^-, (23) 

z w — 1 1 w — i 



where w = uj/ojq. So the change of 11^ depends on uj. 

Let us see what a value AII^ will take when the cycling angular frequency of EWs 
is equal to the angular frequency of GWs, 

OJQ = UJ, (24) 

called the resonant condition. Taking the limit — 1 in Eq.(23) yields a constant 
value 

Atf = (25) 

which has only contribution from the + mode. This is the known result by Cruise [9] . 
Let us discuss other special cases of Eq.(23). 



^7 




Figure 3: The AIl^ as an oscillating function of w when A(q = 0. AII^ has a maximal value 
of 1.864A^ at zu = 1.434 and a minimal value of —0.643^0 at ccj = 0.743. Also notice that: 
(1) AIl^ = At^, for ZJ7 3> 1 and zu equals half integer; (2) AII^ = 0, for zu equals integer; (3) 
An2 = 0, for ui^Q. 
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Figure 4: The AII^ as a function of zu when yl^ = 0. AII^ has a minimum —1.585/1® at 
zo = 1.036. 




(1) If the GWs are given such that ^0 = 0, i.e., there is only the x mode, 

A 2 

Atf = ^ (1 - cos {27rw)) (26) 
2 w — 1 

which is plotted in Fig. 3 as a function of zu. It shows that AII^ can be both positive 
and negative, depending on zu. A maximum value of AII^ is achieved at ~ 1.434. 

(2) If the GWs are given such that = 0, i.e., there is only the + mode, 

Atf = -^sin(27r.^)^-, (27) 

which is shown in Fig. 4. A minimum value of AII^ is achieved at ti7 ~ 1.036. 

(3) If A^ = Aq) = A, where A is real, as is likely the case for relic gravitational 
waves, then one gets 

An^ = ^{w - tx7COs(27rcc7) - sin(27rro)) ^ ^ , (28) 
2 ZU"' — 1 

which is shown in Fig. 5 that is the compound of Fig. 3 and Fig. 4. 
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Figure 5: The AIl^ function of zu when = Aq 
w = 1.546 and a minimum —1.802 at zu = 0.889. 



A. AlP has a maximum 1.842 at 




Figure 6: The relation between AII^ and zu for circularly polarized GWs. The value of AII^ 
will always less than A/ 2. 

Instead of a linearly polarized GWs in Eq.(19), consider the case of circularly 
polarized GWs with = iA^ = A, 



h,^ = A cos (p, hi^ = A sin cp. 
By similar calculations, one has the relevant Christoffel components 

= irAcoscp/Xg, 

= irAsincp/Xg, 



(29) 



Integrating Eq.(21) yields 



, „o Azu sin (Inzu) 
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(30) 
(31) 



which is shown in Fig.6. 

B. GWs travelling along the positive x^-axis 

Different from the above, consider a plane GWs travelling along the positive x^- 
axis. The wave vector is = 27r/Ag(l, 1,0,0). Then, by Eqs.(8) and (17), the phase 
of GW in the torus is 



(f) = -2Trx^/Xg + 2Trx^/Xg = -(27rs + sin (27rs))a;/a;o. 



(32) 
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The metric is now 
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Similar calculations give the relevant ChristofFel components 

r|o = vrA^sin^/Ag, 
rii = - TT sin (j)/\g, 

and the change in 11^ around one circuit of the path 

27r2rA 



(33) 



/ (l + cos(27rs))sinf— (27rs + sin(27rs)))(is. (34) 
Jo \i^o J 



Integrating Eq.(34) gives rise to 



An^ = sin^ (ttzu) 



(35) 



which is oscillates between A^ and 0. Note that the Aq) has no contribution. 

In the case of circularly polarized GWs, Eqs.(33) and (34) should be replaced by 



rio = -'n-Acos(t)/\g, 
r|i = 7rAcos(^/Ag, 



(36) 



and 



An^ 



2Tr^rA 



I (1 + cos (27rs)) cos f—(27rs + sin (27rs))) (is, (37) 



and one has 



Atf = — sin(27rro), 



(38) 



oscillating between A/2 and — A/2. 

c. GWs travelling along the positive x^-axis 
When the plan GWs travel along the positive x^-axis, the metric tensor of spacetime 



is 
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Similar calculations show that the relevant ChristofFel components are 0, and thus 



Mi' 



0. 



(39) 



Thus, the GWs travelling along x^-axis will not change 11^. Therefore, to avoid a null 
result of detection in case of an incident GWs in the x' direction, one may put two 
probes with 90° separation along the annular waveguide. 



Table 1: The case /l^ = 0. The amplitude of (AIl^),, increases linearly with n, as — >• 1. 
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Figure 7: The case of = and n = 10. (An^)^ has a minimal -15.071^® at = 1.00038. 



5 Accumulative effect 

When the EWs pass n cycles along the annular waveguide, the change of may 
be accumulative. This is of practical significance in actual detections. We need only 
consider GWs along the x^— and the x^— directions. 

Firstly, for the linearly polarized incident GWs in the x^— direction, integrating 
Eq.(22) form to n gives the change of 11^ for n cycles 

( An^ )„ = ^ ( 1 - cos (27rnro) ) - sin (27rnro) . (40) 

2 — 1 2 — 1 

In the special case = 0, Fig.7 gives a plot of (An^)„ for n = 10. In contrast with 
Fig.6 for n = 1, (An^)„ is now sharply peaked at ~ 1 with a much larger amplitude, 
as a prominent feature. As given in Table 1, under the resonance condition — > 1, 
the amplitude of (An^)„ increases linearly with n. In fact, this linearly-increasing 
amplitude AII^^-j^ at very large n is also obtained by taking the resonance limit w ^ 1 
of Eq.(40), yielding 

(Atf )„ = --Jl^, (41) 

which is in accord with the result obtained by Cruise [9]. 

The special cases of = and of A^ = A^ are quite similar to each other. 
(An^)„ has, for each given n, both a sharp maximum AII^^j^ at Wmax > 1 and a 
sharp minimum All^j^ at Wmin < 1- As n — ^ oo, the amplitudes AH^^j^, and All^j^ 
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Table 2: The case = 



0. The amplit.Tidcs of cxtrcma of (AIl^),, increase with n. 



n 




A Tt2 / A 

■^min Allj^ijj/A^ 


1 


1.434 1.864 


0.743 -0.643 


10 


1.038 12.027 


0.964 -10.761 


100 


1.0037 114.456 


0.9963 -113.189 


1000 


1.00037 1138.85 


0.99963 -1137.58 


2000 


1.00019 2277.07 


0.99982 -2275.8 


10000 


1.00004 11382.8 


0.999963 -11381.5 



Table 3: The case — A^. The amplitudes of extrema of (An^)„ increase with n. 



n 




^min ^^min/^^ 


1 


1.546 1.842 


0.889 -1.802 


10 


1.055 11.05 


0.982 -20.214 


100 


1.0055 103.504 


0.9981 -205.148 


1000 


1.00055 1028.1 


0.99981 -2051.59 


2000 


1.00027 2055.43 


0.999907 -4109.52 


10000 


1.00005 10274.1 


0.999981 -20549 



increase with n approximately linearly, and their locations Wmax and ^min 

approach 

to 1 from either side, respectively. Fig.8 and Fig.9 give the plots of (An^)„ with n = 10 
for = and for A^ = A(^, respectively. Table 2 and Table 3 list the increase with 
n of the amplitudes of extrema An^„j. and All^j^ for A^ = and for A^ = A^, 
respectively. 

For circularly polarized GWs in the x^— direction, the n— cycle result is 

/ A Atu sin(27rntu) 

= 2(1+^) ■ 

which does not accumulate with n, but vibrates more rapidly than Eq.(31). 

Secondly, for the lineal polarized incident GWs in the x^— direction, the n— cycle 
result is 

{AU\ = A^ sin^ {mrw), (43) 

which has no accumulating effect. For the circularly polarized GWs in the x^— direc- 
tion 

(AUX = (44) 

having no accumulating effect, either. 

So the above detailed analysis on the n-cyclc accumulating effects yields the simple 
conclusion: Only linearly polarized incident GWs in the .x'^-axis has a linearly accu- 
mulating effect of rotation of PV of EWs in the limit — > 1. Theoretically, in order 
to experimentally obtain a maximum effect of n-cycle accumulation, the circling EWs 
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Figure 8: The case of A© = and n — 10. (An^)„ has a maximum An^^^. = 12.027A,g, at 
uj = 1.038, and a minimal An^^„ = -10.761A^ at = 0.964. 
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Figure 9: The case of = A® = A and n = 10. (AII^ 
at ti7 = 1.055, and a minimal All 
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in the waveguide should be running so that n is as large as possible. Of course, due to 
attenuation of EWs in the actual waveguide, for a given waveguide made of conducting 
metal, such as copper, an input beam of EWs in the waveguide can run only a finite 
number of turns around the torus. The maximum value of n is approximately equal 
to the quality factor Q, mainly determined by the conducting metal employed and the 
pump resonances. For instance. Cruise's group [10] has used copper for the waveguide, 
and the measured value of the quality factor Q ~ 2000. As for the selective response 
of the detector to the particular x^— direction of the incident GWs under the resonant 
condition uj ~ wq, it is a problem for gravitational radiations from certain sources, 
since they generally exist for a finite short period of time (from minutes to hours) and 
have some fixed direction of propagation. But for the RGWs as the detection object, 
it is not a problem at all, as they consist of various modes in all directions and of all 
frequencies, moreover, they are a stochastic background existing everywhere and all 
the time. Therefore, the RGWs serve as a natural object of detection. What one needs 
to do is to set up a convenient position of the torus and to fix the the cycling angular 
frequency = v/R of EWs around the waveguide. There are always modes of the 
RGWs with the direction and the angular frequency a; 2± a;o. 



6 Detecting capability for very high frequency 



Let us examine the capability of the waveguide detector built up by Cruise's group [10] , 
particularly in regards to the RGWs. Consider the favorable case of GWs travelling 

along the x^— direction. Since the rotation AIl^ is small, it is equal to the angle a 
rotated, a ~ AII^. This angle can be measured by the electric probe. In general, the 
detector sensitivity will be limited by the thermal noise in the electronic amplifiers. It 
has been found that [10] the minimum detectable angle of rotation 



where / is an efficiency factor of the probe, which transfers electric signals to the 
following electronic amplifiers, k is the Boltzmann's constant, T is the amplifier noise 
temperature, and B is the detector bandwidth in hertz. Thus, for a constant amplitude 
on the time scale ~ Q/i^o, during which the EWs travel Q turns around the loop, by 
Eq.(41), the minimum detectable amplitude hmin of the GWs is 



where the input power Pi,i is related to the circulating power P by Pj„ = P/Q. Here 
Q is the quality factor of the waveguide. For a random signals of GWs with amplitude 
varying considerably over the time scale ~ Q/i^o, the minimum detectable amplitude 
is 



RGW 




(45) 



h. 




(46) 




(47) 



1 4 




5 6 7 H 9 10 II 



Log,, V [Hz] 

Figure 10: The spectrum h{i>,r]H) with different parameter (3^ for (3 = —1.8. 

since the angle a of rotation accumulatively increases as a oc ^/Q as for a random 
walk. 

The waveguide detector is used to monitor the GWs of frequency ~ 10^ Hz, which 
primarily come from the stochastic background of RGWs with a very broad frequency 
range (10^^^ ~ 10^°) Hz [14, 17, 18]. The RGW spectrum h{u,r]H) as given by Eq.(6) 
in a frequency range > lO"^ Hz depends sensitively on the reheating parameter (3s- The 
spectra for three different values of Ps = 0.5, 0, —0.3, respectively, are given in Fig.lO 
for a model with /3 = —1.8, r = 0.22, and 0,a = 0.75. A larger Ps has a lower amplitude 
in the range (10'' ~ 10^) Hz, however, around v > 10^ Hz the spectrum begins to 
increase considerably. Therefore, if the detector is capable of accurately detecting the 
RGWs signals, it will, in principle, be able to constrain the model parameters /3 and 
/3s, and distinguish different models of reheating during the early universe. 

Other sources of GWs, such as binary neutron stars or black holes, merging of 
neutron stars or black holes can produce GWs, but the frequency is much lower than 
10*^ Hz [12]. Thus they are not to be detected by the waveguide detector discussed here. 
There might be other astrophysical processes, which can give rise to high frequency 
GWs [22] . Thermal gravitational radiation of stars can generate GWs at most probable 
frequencies ~ 10^^ Hz. But this frequency is too high for the waveguide detector 
[22] [23]. The predicted graser beams in interstellar plasma can generate GWs with 
"optical" frequencies > 10^^ Hz, still too high for the waveguide detector [24]. The GW 
radiation from primordial black holes of small mass M < 10~^Mq can generate GWs 
with frequencies > 10^'^ Hz, which is too high for the waveguide detector. Besides the 
rate of event is very low ~ 5 x lO^^events/year/galaxy [25]. Therefore, our primary 
object of detection is RGWs, whose spectrum in the current accelerating universe is 
derived in Refs.[14, 17]. 

What the waveguide detector actually detects is the root-mean-square (r.m.s.) am- 
plitude of RGWs per 

HzV2 at a given u, which can be written simply as [18] 



where h{i') denotes the value of the spectrum h{u,r)H) given in Eq.(6). Since the 
waveguide detector works around the frequency 10^ Hz, so we need to examine the 
value of h{u,'qH) around this frequency predicted by our calculation [14, 17]. For a 
cosmological model with the tensor/scalar ratio r = 0.22, the dark energy = 0.75, 
and the reheating 13s = 0.3, one directly reads from Fig.l the values h{i') ~ 10^^^, 10~^^ 
for the values of the inflationary parameter (3 = —1.8, —2.02, respectively. So the 
corresponding r.m.s amplitude per Hz^/"^ at v = lO^Hz is then 



(10-^^ IQ-^^) Hz-^/"^, (49) 



for the two values of (3, respectively. On the other hand, the detector sensitivity can be 
improved by using the cross correlation of two or multiple detectors. From a short run 
of some 4 seconds of the two detectors, Ref.[10] gives the cross correlation sensitivity, 

5 X W-^^Hz-^/\ (50) 

which is within a factor 4 of the predicted sensitivity given that parameters Pin = 
69mW, T = 300K, Q = 2000, {ad)/P = 0.5, and / > 0.9. By comparing the prelim- 
inary experimental result in Eq.(50) with the predicted values in Eq.(49), it is clear 
that the predicted value of RGWs in the model f3 = —1.8 is lower than the prototype 
detector sensitivity by 17 orders. As has been analyzed in Ref.[10], the detector sen- 
sitivity of the current detector could be improved by a factor of 10^ ~ 10^, through 
optimization of the transducers, use of cryogenic amplifiers and multiple detector cor- 
relation. But even with these possible improvements, still that will be some 12 orders 
short to able to measure the predicted amplitude of RGWs in Eq.(49). 

An interesting feature of the spectrum of RGWs is that it has a higher amplitude 
in lower frequencies. This may be suggestive for new ways of enhancing the chance 
of detections. As is seen from Eq.(ll), if one increases the radius R of the annular 
waveguide, say to from 1 meter to 100 meters, the frequency of GWs to be detected 
will subsequently be reduced to a low value ~ 5 x 10^ Hz, at which the spectral 
amplitude h{u,riH) increases by a factor ~ 10^, as is seen from Fig.l. The r.m.s 
amplitude per Hz-^/^ will be h{y)l\fv ~ 10~^^ Hz~-^/^. Now, this is only some 8 orders 
lower than the detector sensitivity of the improved device. Therefore, according to 
our calculation of RGWs, it is unlikely to detect signals of RGWs, using the annular 
waveguide detectors as it stands today. But enlarging the radius i? will enhance the 
detection probability considerably, of course, at the price of a larger sum of cost and 
a more complex construction. By the way, note that LIGO is still unable to detect 
the RGWs by 2 orders of magnitude even it has achieved its design sensitivity [17]. 
Moreover, theoretically, there are possibilities that the waveguide detector can detect 
signals from other kinds of sources of GWS with a much improved sensitivity. 

7 Conclusions 

From the calculations of the rotation of PV of EWs, it is found that the detector only 
essentially responds to the linearly polarized RGWs travelling in the x^-axis under the 



resonant condition. Both the circularly polarized RGWs travelling along any direction 
and the linearly polarized RGWs travelling in the x^- or x^-axis have not observable 
effect. But these propose no problem for the RGWs as the object of detection. 

From our analysis comparing the spectrum of GRWs with the detector sensitivity, 
the RGWs in the accelerating universe have a very low amplitude and are not possible 
to detect using this current detector. The gap between them is some 17 orders of 
magnitude under the current experiment conditions. Even with the improvements 
on the current detector system as planned in [10], there will still be a gap of 12 
orders of magnitude. Examining the detector itself, Eq.(46) and Eq.(47) tell that 
the sensitivity of the detector can be directly improved by several means as follows: 
(1) The use of cryogenic devices at lower temperature T of the environment, i.e., 
to reduce thermal noise of the amplifiers; (2) Increasing the quality factor Q of the 
waveguide, so that the EWs can travel more number of turns around the loop path; 
(3) Enhancing the input power of the EWs into the waveguide; (4) Using multiple 
detectors, whose correlation can improve the sensitivity of the detector. On the other 
hand, the shape of the RGWs spectrum ^(z^) is such that its amplitude is higher in 
lower frequencies. Therefore, it may be more promising to detect the RGWs in the 
relative lower frequency range. For instance, if the radius of torus is increased to 
R = 100 meter, the detecting frequency i^^ ~ 5 x 10^ Hz, and the gap will reduced 
down to 8 orders of magnitude. An overall estimate is that significant improvements 
of the current prototype detector are needed for a possible detection of RGWs by the 
waveguide. 
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